Excitation spectra of transition metal atoms on the Ag (100) surface controlled by 

Hund's exchange 
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We report photoemission experiments revealing the valence electron spectral function of Mn, Fe, 
Co and Ni atoms on the Ag(100) surface. The series of spectra shows splittings of higher energy 
features which decrease with the filling of the 3d shell and a highly non-monotonous evolution of spec- 
tral weight near the Fermi edge. First principles calculations demonstrate that two manifestations 
of Hund's exchange J are responsible for this evolution. First, there is a monotonous reduction of 
the effective exchange splittings with increasing filling of the 3d shell. Second, the amount of charge 
fluctuations and, thus, the weight of quasiparticle peaks at the Fermi level varies non-monotonously 
through this 3d series due to a distinct occupancy dependence of effective charging energies U e fj. 



Strongly correlated multi-orbital quantum systems 
present a classical yet unsolved problem in solid state 
physics which appears in various systems ranging from 
periodic solids to isolated atoms on surfaces. In 
nanoscopic systems, atomic scale control of correlated 
electrons holdspromises for novel modes of information 
processing 0, 0] and offers possibilities to understand 
fundamental quantum effects such as itinerant electron 
magnetism, the competition of local with non-local mag- 
netic interactions |3H5| or the transition from isolated 
atoms to extended solids [sHU . Here, the multi orbital 
nature of realistic nanoscale transition metal (TM) struc- 
tures is generally believed to control physical properties 
such as magnetic anisotropies @, magnetic excitations 
Q or Kondo temperatures [loj . 

While there has been huge the experimental progress 
in fabricating nanoscale correlated electron systems, our 
theoretical understanding of the physical mechanisms de- 
termining their electronic and magnetic properties is still 
at a rather basic qualitative level: Often experiments 
on transition metal impurity systems are interpreted in 
terms of Kondo type models where a generalized spin is 
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coupled to a bath of conduction electrons |M 
in terms of single orbital Anderson models [12h1 
ever, links between these models and realistic systems are 
usually ambiguous, difficult to establish [20| and can typ- 
ically only be made a posteriori It is often largely 
unclear which microscopic degrees of freedom are active 
at a given energy scale and how their contribution in 
excitation spectra measured with different spectroscopy 
techniques can be disentangled. 

In this letter, we consider the series of isolated Mn, Fe, 
Co and Ni adatoms on the Ag (100) surface and explain 
their excitation spectra. We report on valence electron 



photoemission experiments revealing a complex evolution 
of the electronic spectra through this series: we find a 
monotonous decrease in the splitting of higher energy 
features and a non-monotonous variation of low energy 
spectral weight. By means of first-principles calculations 
we explain the photoemission results and show that both 
observations trace back to Hund's exchange. First, the 
splitting between final state multiplets with different spin 
decreases monotonously due to a monotonous reduction 
of effective exchange splittings J(n,f — nj.) with increas- 
ing fillin g of the 3d shell. Second, the effective charging 
energies [MS] U eS (n) = E{n + 1) + E(n - 1) - 2E(n), 
with E(n) being the ground state energy of the impu- 
rity adatom with n electrons, vary due to Hund's ex- 
change in a strongly non-monotonous way from Mn to 
Ni. Therefore, the amount of charge fluctuations and the 
weight of quasiparticle peaks at the Fermi level evolves 
non-monotonously through this 3d series. We find siz- 
able charge fluctuations and mixed valence behavior for 
Fe and Ni. In contrast Mn and Co come closer to a multi- 
orbital Kondo limit with a generalized impurity spin be- 
ing coupled to a bath of conduction electrons and less 
charge fluctuations. 

For our experiments, we prepared the Ag(100) sub- 
strate by the standard procedure and monitored its 
crystalline quality by low energy electron diffraction 
(LEED). The LEED pattern quality was very high, with 
sharp diffraction spots on low background. Isolated TM 
adatoms were obtained by depositing Mn, Fe, Co and 
Ni atoms at a substrate temperature of 20 K (statistical 
growth regime). The TM coverages were calibrated by 
a quartz microbalance and here we conventionally define 
the coverage according to mass equivalent of a nominal 
monolayer. The sample temperature was maintained at 
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Figure 1. (Color online) Valence band spectra of Mn, Fe, Co, 
and Ni adatoms (from top to bottom) of Ag (100). (a) Ex- 
perimental photoemission spectra (b) Theoretical spectra ob- 
tained from QMC calculations at /3 = 20eV~ with imaginary 
time discretization of At — 0.125 via analytical continuation 
[2~il ]. The 3d shell occupancies used in the simulations are 
n = 5.0 for Mn, n = 6.0 for Fe, n = 7.8 for Co, and n = 8.4 
for Ni. 



20 K during the photoemission measurements. The pho- 
toemission spectra were measured with a photon energy 
of 120eV which corresponds to the Cooper minimum of 
the Ag 4c? photoionization cross section. Under these ex- 
perimental conditions, there is an enhancement of the sig- 
nal of the impurity (coverages in the range of 2 x 10~ 2 to 
1CP 3 monolayers) with respect to the one of the host sur- 
face. The photoemission experiments were performed at 
the SuperESCA beamline at the ELETTRA Synchrotron 
Radiation facility, with an overall energy resolution bet- 
ter than 40 meV. 

To explain the experimental spectra we performed 
density functional theory (DFT) calculations using the 
generalized gradient approximation (GGA) [25| as im- 
plemented in the Vienna Ab-initio Simulation Package 
(VASP) (2(| with the projector augmented waves basis 
sets (PAW) [13, H3- In these calculations, single transi- 
tion metal atoms on Ag(100) have been modeled using 
4x4 surface supercells with slab thicknesses of five layers. 
The crystal structures have been relaxed until the forces 
acting on each atom were less than O.OleV/A. This yields 
hybridization functions of the adatoms, which are then 
used to define five orbital Anderson impurity models with 
Coulomb interactions given through the Slater integrals 
F° = U, F 2 = 14/(1 + 0.625) J and F 4 = 0.625.F 2 [H 



with U = 3eV for Mn and Fe (U = 5eV for Co and Ni) 
as well as Hund's exchange J = 1 eV. As the occupan- 
cies n of the 3d impurity orbitals are not exactly known, 
they are kept as free parameters. The impurity models 
were solved usin g th e Hirsch-Fye Quantum Monte Carlo 
method (QMC) [30( (keeping the density-density part of 
the local Coulomb interaction) as well as exact diagonal- 
ization (ED) (3l| , In this way, we obtain the adatom 
spectra including electron correlation effects in QMC as 
well as a detailed insight into the atomic multiplet struc- 
ture via ED. 

Fig- Q] (a) shows the experimental photoelectron en- 
ergy distribution curves in the valence band of isolated 
Mn, Fe, Co and Ni atoms on the Ag(100) surface. The 
curves are difference spectra between the clean Ag surface 
and the surface covered with a few adatoms and thus cor- 
respond to the contribution from 3c? impurity electronic 
states (32] . We observe a remarkable evolution of the im- 
purity spectra through this series of TM adatoms: Mn 
possess one structure (labeled 1) at binding energy (BE) 
of 3.25 eV; Fe has two structures, one (1) at 2.32 cV BE 
and the other one (2) near the Fermi Energy (Ep); Co 
has one (1) broad structure at 2.57 eV BE, one structure 
(2) at 0.8 eV BE, and one structure (3) close to the Ep; 
and Ni has one broad structure (1) at 0.35eV BE. 

We start with reconciling these experimental results in 
the context of a generalized Kondo description: For Mn, 
the spectral peak at — 3.25eV and virtually no quasipar- 
ticle peak at the Fermi level could be well in line with 
Mn acting effectively as a spin S — 5/2 Kondo impu- 
rity. Indeed, this would be very similar to the situa- 
tion found for Mn impurities in bulk Ag, which has been 
derived from photoemission spectroscopy and measure- 
ments of the magnetic susceptibility [5|. The virtually 
absent quasiparticle peak would then be well understand- 
able as the large spin S = 5/2 leads to very low Kondo 
temperatures [lOj. With increasing filling of the 3c? shell 
the impurity spin should be gradually reduced and the 
spectral weight of the quasiparticle peak near the Fermi 
level should be growing exponentially. Indeed, Fe, Co, 
and Ni exhibit spectral weight near the Fermi level but 
the shape and weight of these low energy spectral varies 
very non-monotonically through the series of Fe, Co and 
Ni. In particular, we do not find a monotonous increase 
of the quasiparticle spectral weight as would be expected 
in spin-only Kondo models [Io| . Thus, we conclude that 
additional degrees of freedom must be responsible for the 
experimentally observed evolution of the adatom spec- 
tra. To pinpoint these degrees of freedom and to explain 
the spectra theoretically, we present calculations combin- 
ing density functional theory (DFT) and quantum many 
body methods in the following. 

The DFT calculations show that all transition metal 
adatoms adsorb to high symmetry positions continuing 
the Ag lattice, i.e. sitting in the center of a square of Ag 
atoms. The adsorption height above the surface differs 
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Figure 2. (Color online) Spectral functions for (a) Mn, (b) Fe, 
(c) Co and (d) Ni impurities obtained from QMC calculations 
at P = 20eV _1 with imaginary time resolution of At = 0.125 
via analytical continuation [2J| (thick black lines). The 3d 
shell occupancies used in the simulations are n = 5 for Mn, 
n = 6 for Fe, n = 7.8 for Co, and n = 8.4 for Ni. Additionally 
spectra obtained by ED are shown (orbitally resolved in a) 
and b); total spectra for different filling of the 3d shell in c) 
and d)). 



only little, from Mn at 1.30 A to Ni at approximately 
1.4 A. In line with the similar adsorption geometries the 
hybridization functions, Im A(w), of the adatoms are 
similar for all adatoms 32]. Most importantly, the hy- 
bridization function is rather featureless for all adatoms 
in the energy region between —3 eV and +1 eV. Thus, the 
complex evolution of the spectra observed experimentally 
also cannot be a single particle hybridization effect. 

The spectral functions of Anderson impurity models 
obtained from QMC and ED are shown in Fig. [T] (b) 
and in Fig. [5J In agreement with our experiments, the 
Mn spectrum consists mainly of one peak far below the 
Fermi level for 3c? shell fillings n « 5. In fact, already a 
diagonalization of the atom in the crystal field (A c /) of 
the surface shows the basic structure found in the exper- 
iment (Fig. [5^). Thus, a low energy description of Mn 
on Ag (100) in terms of a spin S = 5/2 Kondo model is 
well in line with our results. This is indeed similar to the 
case Mn in bulk Ag [5j and also in agreement with DFT 
calculations for Mn on Ag(100) [33, 34 1. Comparison of 
the QMC calculations to the experimental spectra reveals 
good agreement also for the Fe, Co, and Ni adatoms. 
We thus use the QMC results to understand the physical 
mechanisms behind the evolution of the spectra in the 
series of 3d adatoms. 

Multiplet splittings and Hund's exchange. For Fe with 
n = 6, our calculations show a broad peak around — 3eV 
and a relatively narrow peak right below the Fermi level, 



Figure 3. (Color online) Orbitally resolved spectral functions 
of Co (n = 7.8 (a) and n = 7.9 (c)) and Fe (n = 6.0 (b) and 
n = 6.4 (d)) impurities as obtained from QMC. 



which reproduce the experimental features 1 and 2, re- 
spectively. The broad satellite around —3 eV appears in 
all orbitals and is found also in the ED calculations (Fig. 
EJd) . It can be identified as d 6 — > d 5 ionization peak. An- 
alyzing the orbitally resolved spectral function shown in 
Fig. [3jb) we see that the experimental feature 2 stems 
from the d x 2_ y 2 orbital. This orbital has the occupancy 
0.8 in contrast to the other orbitals which are approxi- 
mately half filled. Therefore, the feature 2 cannot be a 
quasiparticle peak due to spin-only Kondo physics. In- 
deed, this peak appears already in an ED description of 
Fe in the crystal field (Fig. [2] (b)) of the surface, which 
shows that it corresponds to a d 6 — > d 5 transition. The 
spectral features 1 and 2 of Fe, thus, stem both from ion- 
ization processes of the impurity. Our ED calculations 
further show that the energy separation of these peaks 
traces back to different d 5 final state multiplets: S = 5/2, 
L = for feature 2 near Ep and higher energy multiplets 
such as S = 3/2, L > for feature 1. The splitting be- 
tween these multiplets can be understood as an effective 
exchange splitting ~ J(ji| — n^) due to Hund's exchange. 

The experimental Co spectrum consisting of three 
peaks is well reproduced in our QMC simulations for 
n = 7.8 (Fig. [2t). Analyzing the orbitally resolved spec- 
tral function shown in Fig. [3] (a, c) and the corresponding 
occupation matrices we find that the d XZiVZ and d x i_ y i 
orbitals are almost fully occupied and mostly responsible 
for the peak at — leV, which corresponds to the feature 
2 in the experimental spectra. Again, the spectral weight 
further below the Fermi level (feature 1) traces also back 
to an ionization process, here d 8 — > d 7 , with higher en- 
ergy final state multiplets. As in the case of Fe, the split- 
ting between the features 1 and 2 can thus be understood 
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as atomic multiplet effect due to Hund's exchange. For 
Co, the effective exchange splitting reduces as compared 
to the case of Fe by an amount on the order of J. This 
explains why the separation of the peaks 1 and 2 is 0.5 
eV smaller for Co than for Fc. 

The spectrum of Ni turns out to consist mainly of 
a broad peak below the Fermi level without clearly re- 
solvable multiplet features. This is qualitatively in line 
with an even further reduced effective exchange splitting 
- J(n t - n±) for Ni. Indeed, in our GGA and GGA+U 
calculations, the Ni adatoms turn out to be nonmagnetic 
on this surface which is in agreement with calculations 
in Ref. [35| and experiments for Ni on Au [3(3]. 

Effective charging energies and valence fluctuations. It 
remains to be explained why feature 2 in the Fe spectrum 
and the whole Ni spectral peak are rather close to Ep 
and what the nature of the peak 3 near Ep in the Co 
spectrum is. Therefore, we discuss the issues of valence 
fluctuations as well as Kondo physics in the following. 

For a fully rotationally invariant Coulomb vertex there 
is a pronounced occupancy dependence of the effective 
charging energies [H, |37} 



U + iJ for n = 5 
U eS (n) » \ U - (3/2) J for n = 6 and n = 9 
U - (1/2) J for n = 7 and n = 8. 



(1) 



Mn has n — 5 and thus the highest U e g which further 
corroborates our conclusion of Mn resembling an atomic 
spin S = 5/2 with nearly frozen valence and is in line 
with the discussion of Kanamori type Coulomb inter- 
actions in Ref. 2l|, |22j. In contrast to the Kanamori 



model, here U e s also varies between the non-half filled 
cases (n = 6, ...,9). Most importantly, we find that the 
d 6 and d 9 atomic configurations yield the smallest U e g 
and are most susceptible to valence fluctuations. This 
gives a hint towards mixed valence behavior of the Fe 
and Ni adatoms, which is substantiated by our QMC re- 
sults. 

For the Fe d\ orbital, there are no well defined upper 
Hubbard bands but only spectral peaks above Ep which 
extend to or even below the Ep (Fig[3]b) and d)). We 
further find that this overall structure of the spectra re- 
mains stable also at larger fillings, like n = 6.4 shown 
in Fig|3][d). In this entire range (6 < n < 6.4) of oc- 
cupancies [38j], the Fe adatoms are in a mixed valence 
situation. 

The experimental spectrum of Ni consists mainly of a 
broad peak below the Fermi level and we find good agree- 
ment between the QMC calculations and the measured 
3d spectra of Ni adatoms for occupancies 8.3 > n > 8.7. 
[39] ] In this range, the Ni spectra obtained from our QMC 
simulations and the experiments are qualitatively more 
similar to the " non- interacting" GGA density of states 
than to the ED spectra shown in Fig|2]d). There are no 
well defined upper Hubbard bands in any of the Ni or- 
bitals but only broad spectral weight distributions above 



Ep which extend below Ep. This points towards a mixed 
valence situation for also for Ni. 

In this respect Fe and Ni are very different from Co: 
For Co, our experiments (feature 3) and calculations 
show a quasi particle peak at the Fermi level, which is 
well separated from clearly formed upper and lower Hub- 
bard bands. There are, thus, less charge fluctuations for 
Co on Ag (100) and so this system comes closer to the 
(multiorbital) Kondo limit. Therefore, our results con- 
firm the interpretation of low energy resonances in STM 
spectroscopy experiments of Co on Ag(100) in terms of 
a Kondo effect 40] . The QMC results further show that 
all Co orbitals are involved in the quasiparticle resonance 
(Fig. [3]a,c). Thus, excitations of the orbital degree of 
freedom must be available at energies on the order of our 
simulation temperature 1//3 = 0.05 eV. This is similar to 
the case of Co on Cu (111) 

For Co we find best agreement of calculated and mea- 
sured spectra in the range n — 7.8 — 7.9 (Figj3] a,c). 
Thus, the Co is closer to a c? 8 than to a d 7 configura- 
tion which supports recent coupled cluster calculations 
[I2I ]. We note that a prediction of the Co valency based 
on DFT type approaches can be misleading: LDA+U 
calculations with 2 eV< U < 5 eV yield an occupancy of 
the 3c? shell between n = 7.0 and 7.2. 

The spectral function of Mn (a group VH-element) 
could be well understood assuming a filling of n « 5 for 
the Mn 3c? orbitals. However, noble metals like Cu (group 
Xl-elements having 4 electrons per atom more than the 
corresponding group VII elements) have an almost full d 
shell (i.e. n « 10) due to one electron from the 4s orbitals 
being promoted to the 3c? orbitals. If this promotion of 
one electron from the 4s to the 3c? orbitals would occur 
homogeneously, the 3c? occupancy should increase by 1.25 
electrons between each two atoms of the 3c? series under 
investigation. Such an increase is in line with our results 
and the mixed valence behavior for Fe (6 < n < 6.4) and 
Ni (8.3 > n > 8.7) but not with Co which comes closer 
to the Kondo limit of nearly frozen c! 8 valence. 

In summary, our joint experimental and theoretical 
study shows that Hund's exchange controls the physics 
of 3c! adatoms on the surfaces of Ag (100). It fosters the 
formation of multiplcts, determines multiplet splittings 
and modulates effective charging energies. Our results 
show that any realistic description of magnetic nanosys- 
tems should account for these manifestations of Hund's 
exchange. Particularly the mixed valence behavior of Fe 
and Ni challenges discussions of transition metal based 
nanomagnetic structures or impurity systems in terms of 
spin-only models. The situation is complex and challeng- 
ing. Further spectroscopy studies, including photoclcc- 
tron spectroscopy with higher energy resolution and as a 
function of temperature, will be useful to fully describe 
the nature of general excitation spectra of magnetic im- 
purities on surfaces. 
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